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ABSTRACT 

We continue to see a range of values for the Hubble constant obtained from gravitationally 
lensed multiple image time delays when assuming an isothermal lens despite a robust value 
from the HST key project (72 ± 8 kms -1 Mpc -1 ). One explanation is that there is a vari- 
ation in Hubble constant values due to a fundamental heterogeneity in lens galaxies present 
in groups, i.e. central galaxies with a high dark matter surface density, and satellite galaxies 
with a possibly stripped halo, low dark matter surface density, and a more central concen- 
trated matter distribution. Our goal is to see if a variety of group interactions between the 
most massive group members can result in significant changes in the galaxy density profiles 
over the scale probed by strong lensing (<15 kpc). While stripping of the outer parts of the 
halo can be expected, the impact on inner regions where the luminous component is impor- 
tant is less clear in the context of lensing, though still crucial, as a steepened density profile 
within this inner region allows these lens systems to be consistent with current HSTAVMAP 
estimates on Hq. We employ the particle-mesh code SUPERBOX to carry out the group in- 
teraction simulations. An important advantage of using such a code is that it implements a 
fast, low-storage FFT-algorithm allowing simulations with millions of particles on desk-top 
machines. We simulate interactions between group members, comparing the density profile 
for the satellite before and after interaction for the mass range of 10 11 to 10 13 Mq. Our inves- 
tigations show a significant steepening of the density profile in the region of ~ 5-20 kpc, i.e. 
that which dominates strong lensing in lens galaxies. This effect is independent of the initial 
mass-to-light ratio. Additionally, the steepening in the inner region is transient in nature, with 
consecutive interactions returning the profile to an isothermal state within a timeframe of ~ 
0.5 - 2.0 Gyr. This factor may help explain why lens galaxies that produce lower values of Hq 
(i.e. those with possibly steeper profiles) are far fewer in number than those which agree with 
both the HST key project value for Ho and isothermality, since one would have to observe the 
lens galaxy during this transient steepened phase. 

Key words: Methods: N-body simulations - galaxies: structure - galaxies: interactions - 
galaxies: haloes - Gravitational lensing 



1 INTRODUCTION 

A range of values for the Hubble constant is obtained from 
different strong gravitational lens systems with measured image 
time dela ys, assuming isothermal models for the lensing galax- 
ies (e .g. iBiggs et all 1 19991 : iGil-Merino. Wisotzki & Wambsganssl 
2002). The estimates are typically lower than the values from 
the HST key project and WMAP which are consistent with each 
other at 72 ± 8 kms" 1 Mpc" 1 and 72 ± 5 hns" 1 Mpc~\ re- 
spectively i Free dman et al.i 2001 : Sp ergel et al.ll2003h . Since this 
Hubble constant value is considered robust, this variation raises 
questions as to the relative extent of the galaxy halo and the lu- 
minous component since the derived values fundamentally de- 



pend on the assumed lens potential and matter surface densi- 
ties. A number of differently motivated studies have found that 
the inner regions of lens galaxies exhibit is othermal density pro- 
files (e.g. iRusin. Kochanek & Ke eton 2003 ), including rec ent re- 
sults from the Sloan SLACS survey iKo opmans et al.ll2 006) which 
found t] = 2.01 to'03 (lc; ptot oc r~ v ) when studying 15 early- 
type galaxies. Moreover, studies inferring the Hubble constant from 
10 time-delay lens galaxies produced values in agreement with 
the HSTAVMAP value, while still being consistent with simu- 
lated tim e-delays from N- body galaxies which were isothermal 
in nature dSaha et alj|2006h . Even a relatively simple analysis us- 
ing the existing time delay system s indicates a consistency with 
isothermality JDobke & King|l2006h . Stellar dynamical studies also 
point towards isothermal density profiles in the inner regions of 
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galaxies (e.g. Rix et al. 1997 ;_ Romanowskv & KochanelJ 1 19991 ; 

2d). 



iGerhard et al J200 ll ; iTreu & Koopmansll2002; 



Despite this apparent trend towards isothermality, we still ob- 
serve a number of systems where it is hard to reconcile Ho val- 
ues with the HST key project value when assuming an isother- 
mal model. Modelling of the systems PG1 115+080, B 1600+434, 
HE2149-2745, and SBS1520+530 has shown that these systems 
prefer more centrally concentrated density profiles and hence non - 
isothermal models dKochane k 2002; K ochanek & Sche chter 2004). 
In the specific case of PG1 115+080, detailed modelling done by 
combining lensing, stellar kinematic and mass-to-light ratio con- 
straints, in order to build a two-component model, shows a mass 
density profile significantly steeper than isothermal at 77 = 2.35 ± 
0.15, where p to t oc r~ v . However, the derived value for the Hubble 
constant came to Ho = 59 ly 2 kms -1 Mpc -1 , implying the need 
for an even steeper profile to agree with the HST key project value 
dTreu & Koopmans 20 02rJ). 

iKeeton & Zabludofl d2004l) have shown that lens models of 
four-image systems overestimate the Hubble parameter by up to 
15% when neglecting the effect of the contribution from other 
group members on the potential of the lens galaxy, with an even 
higher discrepancy for two-image systems. Hence accounting for 
environment would further reduce the derived value of Ho for a 
fixed primary lens model, requiring the slope of the primary lens 
to be further steepened to obtain consistency with the HST key 
project value. Additionally, outside the group environment, line- 
of sight structures can provide contributions to the lens potential 
(Mo mcheva et alj2006h . 

While isothermal systems certainly appear to be the more 
prevalent, an explanation of this sub-sample of lens galaxies re- 
quiring non-isothermal profiles is not clear. One suggestion fo- 
cuses on the possible role that the group may have on the lens 
gala xy, specifically the eff ect of tidal stripping of satellite galax- 
ies ( Koch anek etal] 120061) as recently observe d in clusters us- 
ing galaxy-galaxy lensing dLimousin et alj|2006l) . Although a full 
census of lens galaxy environments is n o t available, a number 
of observations dTonrv & Kochanek! 1 19991 ; iFassnacht etal] |2002| ; 
IWilliams et all [2006) seem to show that typical early-type lens 
ga laxies should be membe rs of a group. In the halo model described 
bv lCoorav & Shethl d2002h . one galaxy typically lies at the centre of 
the group halo, while other galaxies are smaller satellites orbiting 
in that halo. In this scenario we have two differing cases for galaxy 
structure. In the first we see the central galaxies with a high dark 
matter surface density and the luminous component making only a 
small fraction of the overall matter content. In the second case, we 
see a satellite galaxy with a possibly stripped halo, and more cen- 
trally condensed core. In essence, it suggests a fundamental hetero- 
geneity in early-type galaxies present in groups, and when applied 
to lens galaxies the distinction between the two surface mass den- 
sities could give rise to a difference in the derived Hubble constant 
were a standard isothermal model assumed. 

This description of the heterogeneity in groups begs a simple 
question; how much can the group members actually effect each 
other in terms of changing their mass and density profiles? In such 
groups, unlike large clusters, there is not so much a single dom- 
inant central galaxy, but rather a power sharing scenario between 
two or more approximately equivalent galaxies. While interactions 
between members with large differences in their respective masses 
should be expected to produce large dynamical and structural ef- 
fects in the less massive member (e.g. ultra compact dwarf (UCD) 
interactions; see Fellh auer & KroupJ r2006). the effect of interac- 
tions between more similarly matched group members is perhaps 



not so clear when considering the effect in the context of lens galax- 
ies. Indeed, although c ertain past studies (e.g. lHavashi et afll2003l ; 
Kazantzidis et al. 2004) have investigated the effect of satellite- 
halo interactions, the focus has been on the evolution of substruc- 
ture or of the very inner cusp rather than the effect on cosmological 
parameter derivation from lens galaxies. 

The question posed is a complex one, and the parameter space 
of group interactions is vast. In this paper our aim is to establish 
whether, over the scale probed by strong lensing (<15 kpc), signif- 
icant changes in the density profile of a galaxy on a satellite orbit 
can be induced by interactions with a group galaxy at a central lo- 
cale. While stripping of the outer parts of the satellite's halo could 
be easily expected, the impact on inner regions where the lumi- 
nous component is important is less clear in the context of lensing, 
though crucial, as it is with a steepened density profile within this 
inner region that the time delays for these lens systems can agree 
with current HST/WMAP estimates on Hq. In this way we ascer- 
tain if another mechanism is required to alter density profiles, pos- 
sibly during the group's formation. 

The outline of the paper is as follows; in §2 we briefly intro- 
duce the N-body code used to simulate the interactions, with details 
of the exact parameters of the simulated objects. In §3, we go on 
to present results and discuss their analysis in the framework of 
the problem. Finally, §4 draws some conclusions based upon our 
findings. 



2 THEORY AND SIMULATION SETUP 

2.1 The particle-mesh code SUPERBOX 

Our analysis employs the particle-mesh code SUPERBOX to carry 
out the simulations. The code utilises high resolution sub-grids and 
a nearest grid point (NGP) force-calculation scheme based on the 
second derivatives of the input potential. An important advantage 
of using such a code is the fact that it implements a fast, low-storage 
FFT-algorithm giving the possibility to work with millions of par- 
ticles on desk-top machines or small clusters. 

For each object used in the simulation, there are 5 grids with 

3 different resolutions. For the purposes of our simulation setup, 
we define the size of the inner, middle and outer grid sizes. The 
inner grid is the highest-resolution grid which resolves the centre 
of the galaxy. The middle grid has an intermediate resolution to 
resolve the galaxy as a whole, while the outer grid is the size of the 
whole simulation area (i.e. 'the local universe'), and has the lowest 
resolution. The inner and middle grid remain focused on the object 
throughout the simulation. Each grid has a predefined number of 
mesh-points. In our simulations we use 2 6 mesh-points. Doubling 
the number to 2 7 showed no significant effect on the overall results 
discussed in §3, indicating a high enough level of resolution. For 
a full description of the code and grid structure see lFellhauer et al] 
d200d) . 



2.2 Simulation Setup 

We focus the main body of our investigations on approximately 
imitating the dynamics and parameters of lens systems with a 
known history of lens models with low dark matter surface densi- 
ties, and hence possible steeper density profi les, e.g. PG1 1 15 +080, 
B 1600+434, HE2149-2745, SBS 1520+530 dKochanek! I2OO2I) . and 
as such are clear candidates for possible satellite stripping. We use 
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Figure 1. Left: Contour plot of an NFW halo after undergoing isolated integration. The contours have (logarithmic) magnitude spacing with the innermost 
and outermost contours corresponding to 16 and 36 mag arcsec -2 respectively. The object shows clearly defined symmetrical structure and is ready to be 
used in a simulation scenario. Right: The evolution to dynamic equilibrium as shown by the Lagrangian-radii of the mass shells of 10% to 90% of particles 
in a typical simulation. The flattening of the lines with time shows that an equilibrium state has been reached, and that the object is ready to be used in an 
interaction scenario. 
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Figure 2. Left: Example of a path taken in the x-y plane during a simulation run. The orbit starting position in the x-y plane is (0,-300) kpc, with a tangential 
velocity Vtan = +100 kms— 1 in the x-dimension. This particular run shows the interaction between a 10 11 Mq object after undergoing three passes of a 
central 10 12 Mq object. The run time is 7.3 X 10 9 years. Also shown are the 100 kpc and 150 kpc radial distances from the central potential, beyond which 
the satellite spends 78% and 64% of its orbit time respectively. Right: Typical contour plot for the same 10 11 Mq halo after undergoing two interactions with 
a 10 12 Mq central object. Note the formation of tidal-tail structures in the outer regions. 



the term satellite here to refer to a galaxy undergoing orbital inter- 
actions with another galaxy, whether that be of equivalent or higher 
mass, and not in the sense of e.g. a dwarf elliptical satellite. Some of 
the p arameters used to se t up the model galaxies were tak en in part 
from Imp ev et alj dl998h and lKeeton & Zabludof J J2004) . Param- 
eters such as mass and scale are obviously often model dependent 
and as such were only taken as approximate guide values and sanity 
checks. In actuality, variations of up to an order of magnitude in the 
masses of the model objects had little effect on the overall results; 
more important were the relative masses between the interacting 
bodies. 

For the dark matter halo component of each of the group 
galaxies we use a Navarro, Frenk & White (NFW) profile 
dNavarro."Fr enk & White 1996) while for th e luminous component 
we use a Hernquist profile (Hernquist 1990). Hence, the combined 
density profile for each of the modelled galaxies has the form, 



ptot{r) = piv(r) + pn{r) 

_ Po,n r s>N ^ M H r B , H 
r[l + (r/r s ,jv)] 2 2ty r(r atH + r) 3 

We have the option in SUPERBOX to either define the charac- 
teristic density for the NFW halo, po,N, or provide the total mass, 
Mjv, in which case the code will calculate the density in equation 
CD by means of the following, 

Po,n = ( In(r a w + r) - Infrvjv) ) (2) 

47rr s,jv V r s>JV + r J 

In order to fully define the objects used in the simulations we 
must provide values for the following parameters; NFW total mass 
(Mjv) and scale length (r s ,jv), the Hernquist total mass (Mjj) and 
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Table 1. Table of our initial model parameters for three particular interactions. The first two rows correspond to the object parameters that were used in the 
simulations of Fig.[3]and Fig.f4] We include NFW and Hernquist mass with corresponding scale and cut-off radii, the characteristic crossing time (T cr ) and the 
grid resolutions used in each case. We use the virial radius as an a pproximate guide fo r th e cut-off radii. T he middle grid size is then chosen to encompass this 
value. Parameters for the 10 13 Mq object were taken in part from lMcLaughliij 09991) and lVesperini et alj 420031) 





Initial Object Parameters 




Correspondin; 


g Analytic Potential Parameters 


Grid Resolution 


Massjv.H (M ) 


r a ,Nlr a ,H (kpc) 


r c ,N/r Ct H (kpc) 


T cr (Myr) 


Massjv,// (M Q ) 


r s ,N/r 3 ,H (kpc) 


r c ,N/r CtH (kpc) 


In./Mid./Out. (kpc) 


lO^/lO 10 


10/10 


120/120 


292.36 


io^/io 11 


20/10 


250/120 


50/150/1000 


10 12 /10 n 


20/10 


250/120 


250.58 


10 12 /10 n 


20/10 


250/120 


50/300/1000 


lO^/lO 11 


20/10 


250/120 


250.58 


10 13 /10 12 


560/5 


150/150 


50/300/1000 



scale length (?%,,#)> and finally the cut-off radii for both, r c ,N, r C: H- 
These last two values are used by the simulation during setup to de- 
fine some physical limit to the objects, since both potentials stretch 
to infinity. In general, the virial radius of the object is used as an ap- 
proximate guide for the cut-off radius. In addition, it is necessary 
to define the number of particles that each component (Halo + Lu- 
minous) is modelled with. In all our investigations, both the NFW 
and Hernquist components were modelled using 10 6 particles each. 
We model the central galaxy as a fixed analytic potential (which is 
in effect equivalent to observing the interaction in the rest frame of 
the central galaxy) and observe the effect that the interactions have 
on the satellite. Although the use of a static potential eliminates the 
effect of dynamical friction, in this analysis we are concerned with 
the effect of rapid 'fly-by' interactions and not long term evolution. 
In addition, dynamical friction is more of a factor for many-orbit 
interactions, and we restrict ourselves to simulations lasting only 
a few interactions (i.e. a few pericenters), which is often less than 
two complete elliptical orbits, e.g. Fig.[2](left). Moreover, replacing 
the static potential with a particle based object confirmed that the 
orbits are changed by dynamical friction by < 5 % in these types of 
interaction. The parameters used for the analytic potential are the 
same as for the particle-object and as such are simulated with the 
same expressions, Q} and {2}, above. 

The grid sizes discussed in §2.1 are chosen to compliment the 
dimensions of the object. For the inner grid it is necessary that the 
dimensions cover at least a couple of scale radii. The middle grid 
is chosen to encompass the object, and is set as the virial radius 
of the object plus some additional margin. The outer grid of the 
simulation is always chosen to allow the objects full range of dy- 
namical motion. For this investigation the outer grid, and hence the 
simulation area, was set as IMpc x IMpc. 

Models are set up in virial equilibrium making use of the Jeans 
formula as follows, 



(Vr)i = 



pi(r) 



dr' 



pi(r')dr' 



(3) 



where 4>totai(r') is the combined potential of the NFW and Hern- 
quist components. Having (vf.) we generate velocities with a Gaus- 
sian distribution. 

The objec ts ar e constructed using th e procedure described in 
iKroupal d 19971) and lFellhauer et al.1 {2000). Once the grid sizes are 
chosen, and the parameters of the object are set, we allow the ob- 
ject to be integrated in isolation until relaxation with the grid sizes 
has occurred and the object has reached equilibrium. This is nec- 
essary for all numerical realisations of stationary self-gravitating 
systems because discreteness of the model always implies an ini- 
tial mismatch between kinetic and potential energy. Equilibrium is 
determined by observing the fluctuation of the Lagrangian-radii of 
the mass shells of 10% to 90% of particles present in the simula- 



tion. An example of this is shown in Fig.Q](right). A stable flat line 
shows that no energy is entering or being bled from the system, and 
that the object has reached a stable equilibrium. Once the objects 
have reached this point, they can then be used in the desired simu- 
lation scenario. For example, Fig.[T](left) displays a typical contour 
density plot for a stable NFW that has undergone relaxation and is 
ready to be used in a simulation. Also shown in Fig.|2] (right) is an 
example of the same contour plot after undergoing an interaction 
scenario such as will be detailed in §3. We assume the interactions 
to be spherically symmetric and hence model the interactions in 
the x-y plane of the simulated area, and no prescription of adiabatic 
contraction is included when setting up the objects. 



3 ANALYSIS AND RESULTS 

We have performed a number of simulated interactions covering 
a variety of galaxy masses and initial orbits. The mass range in- 
vestigated is ~ 10 11 to 10 13 Mq. In all interactions the central 
galaxy was located at (0,0) kpc in the x-y plane. We vary the orbit 
of the satellite galaxy by changing both the initial starting position 
and tangential velocity (within the bounds of typical velocity dis- 
persion for group members). For the determination of our results 
we analyse the satellite features after modelling through a signifi- 
cant fraction of the Hubble time. In most cases this involves several 
passes of the central galaxy. An example of a typical orbital interac- 
tion path is displayed in Fig[2](left). The central galaxy lies at (0,0) 
kpc, and the satellite galaxy begins at (0,-300) kpc with a tangen- 
tial velocity of + 100 km s _1 in the x-dimension. We also show 100 
kpc and 150 kpc radial distances from the central potential, beyond 
which the satellite spends 78% and 64% of its orbit time respec- 
tively. This provides some quantitative estimate and understanding 
for the probability of finding such a satellite galaxy away from the 
central galaxy and group centre. 

We highlight results from two particular simulations. The find- 
ings are non-unique to the particular set up and the results were 
found to be quite general for a variety of initial conditions. The 
main results are displayed in Fig.|3]and Fig.|4] We display the den- 
sity profile of the satellite before and after undergoing the inter- 
action, then focus on the region <15kpc as this dominates strong 
lensing. We fit the combined density profiles with a simple power- 
law model of the form, 



ptot{r) = p r~ 



(4) 



where the free parameters are po (some normalisation constant) and 
77, the steepness of the profile. For the fitting procedure itself we use 
equal sized logarithmic bins and employ a nonlinear least-squares 
(NLLS) Marquardt-Levenberg algorithm. 

In Fig.[3] (left) we show the initial combined density profile 
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Initial density profile 
Density profile after 2nd interaction 
Density profile after 3rd interaction 




Initial density profile 
Top solid line: Power-law fit t\ = 2.01 
Density profile after 2nd interaction 
Bottom solid line: Power-law fit t\ = 3.05 




Initial density profile ■ 
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Density profile after 3rd interaction ■ 
Bottom solid line: Power-law fit ti = 2.04 ■ 
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Figure 3. Left: Combined density profile (NFW + Hernquist) before and after a 2nd and 3rd interaction between a 10 11 Mq galaxy and 10 12 Mq object 
modelled with an analytic potential. Initial x-y position and velocity of the satellite were (0,-300) kpc and (+60,0) kms" 1 respectively. Middle: Enlarged plot 
of the density profile in the region <15kpc after the 2nd interaction. Dashed line is the combined density profile before interaction, while dotted line is the 
combined density profile after interaction. Top solid line shows a power-law fit with r\ = 2.01 ± 0.03 i.e. isothermal, while bottom solid line shows a power-law 
fit with rj = 3.05 ± 0.07. The profile has steepened significantly. Right: Enlarged plot of the density profile in the region <15kpc after the 3rd interaction, ~ 
2 X 10 years after the 2nd interaction. Dashed line is the combined density profile before interaction, while dotted line is the combined density profile after 
interaction. Top solid lines show a power-law fit with r\ = 2.01 ± 0.03, with the bottom at r\ = 2.04 ± 0.06. The profile has returned to an isothermal state. 



Initial density profile 
Density profile after 3rd interaction 
Density profile after 4th interaction 




Initial density profile 
Top solid line: Power-law fit t\ = 2.00 
Density profile after 3rd interaction 
Bottom solid line: Power-law fit t\ - 2.70 
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Figure 4. Left: Combined density profile (NFW + Hernquist) before and after a 3rd and 4th interaction between a 10 12 Mq galaxy and 10 12 Mq object 
modelled with an analytic potential. Initial x-y position and velocity of the satellite were (0,-300) kpc and (+100,0) kms" 1 respectively. Middle: Enlarged 
plot of the density profile in the region <15kpc after the 3rd interaction. Dashed line is the combined density profile before interaction, while dotted line is the 
combined density profile after interaction. Top solid line shows a power-law fit with r\ = 2.00 ± 0.02 i.e. isothermal, while bottom solid line shows a power-law 
fit with Tj = 2.70 ± 0.06. The profile has steepened less so than in the case shown in Fig.[3] although the effect is still significant. Right: Enlarged plot of the 
density profile in the region <15kpc after the 4th interaction, ~ 0.5 X 10 9 years after the 3rd interaction. Dashed line is the combined density profile before 
interaction, while dotted line is the combined density profile after interaction. Top solid lines show a power-law fit with r\ = 2.00 ± 0.02, with the bottom at r\ 
= 2.01 ± 0.20. The profile has returned to an approximately isothermal state. 



(NFW + Hernquist) for the simulated satellite galaxy (topmost 
line), and the density profile after the 2nd and 3rd interactions im- 
mediately below. The masses of the interacting bodies were 10 11 
and 10 12 Mq for the satellite and central galaxy respectively. The 
initial mass-to-light ratio of both galaxies was set to 10:1, while 
the starting position was set at (0,-300) kpc, with an initial tangen- 
tial velocity of +60 km s" 1 in the x-dimension. For the purposes of 
this study we are interested in how the interactions in the outer part 
of the halo effect both the luminous and dark matter components 
in the inner part of the galaxy, particularly the region <15kpc as 
this scale encompasses the strong lensing regime. In Fig.[3] (mid- 
dle), we enlarge this region, just comparing the initial and post-2nd 
interaction density profiles. Here, the dashed line is the profile be- 
fore interaction, while the dotted line is the profile after interaction. 



The top solid line shows a power-law fit with r) = 2.01 ± 0.03, i.e. 
isothermal, while the bottom solid line shows a fit with r\ = 3.05 ± 
0.07. We clearly see significant steepening occurring after the 2nd 
interaction over the region of ~ 5 - 20 kpc. The right plot of Fig.|3] 
focuses on comparing the initial and post-3rd interaction density 
profiles. The time elapsed between the two consecutive interactions 
is ~ 2 x 10 9 years, during which the slope of the profile remains 
approximately constant at r\ ~ 3.0. We now see that after the 3rd 
interaction, the slope has returned once more to an isothermal state 
at 77 = 2.04 ± 0.06, and hence the steepening observed after the 2nd 
interaction appears transient in nature. 

We see a similar effect taking place in Fig.[4] except in this 
case we are investigating the interactions between two galaxies of 
similar mass, both being 10 12 Mq. We see that the reduction in the 
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Figure 5. Left: The time evolution of kinetic energy of the bound remnants of the interacting satellite of Fig.[5] Note the three fluctuations in energy resultant 
from the three interactions and that this energy change becomes progressively smaller with time and hence the steepening of the profiles becomes less severe 
after each passage. Middle: The time evolution of potential energy of the bound remnants of the interacting satellite of Fig. [5] Again we see the three energy 
fluctuations, each becomes less severe after each passage. This results in the profile being ultimately unchanged despite a transient steepening during the 
first interactions. Right: The time evolution of potential energy just within the region of 5-20kpc, for the case of the interacting satellite of Fig.[3] We see the 
potential well increase in depth after the interactions, confirming the steepening of the profile, which again lessens after each interaction. 



relative mass necessitates a greater number of interaction events to 
initiate the steepening when compared with the previous case. In 
the middle plot of Fig.[4]the steepened profile after the 3rd inter- 
action can now be fitted with a slope of rj = 2.70 ± 0.06 over the 
range of ~ 6 - 20 kpc. Although less than in the 10 11 vs. 10 12 
Mq case, this is still a significant effect. Once again, after the next 
interaction we observe the profile returning to an approximately 
isothermal state with rj = 2.01 ± 0.20 (Fig.[4] right). This time the 
transient steepening lasts for ~ 0.5 x 10 9 years. 

We now consider the steepening effect seen in Fig.[3]and Fig.|4] 
in terms of energy transfer during the interactions. Displayed in 
Fig.[5]is the time evolution of the kinetic and potential energies of 
the bound remnants for the interactions shown in Fig.[3] During the 
fly-bys of the satellite galaxy, energy gets transformed from orbital 
to internal energy. Particles in the satellite gain kinetic energy and 
those that are loosely bound, mainly in the outer parts, get stripped 
resulting in the steepening of the outer profile. Closer towards 
the inner portion of the galaxy tidal shocks are more important 
than s tripping (e.g. lRead et al.l2006l : lGnedm. Hernquist & OstrikerJ 
1999). Close pericentre passes can generate matter compressions 
possibly resulting in the observed transient fluctuations in the den- 
sity profile. As the mass gets redistributed and the object returns 
once again to equilibrium, the original slope returns to its approxi- 
mate initial state. Notice in Fig.[5]that this energy change becomes 
progressively smaller with time and hence the steepening of the 
profiles becomes less and less severe after each passage, resulting 
in a return to isothermality. With greater time we see no change in 
the profile, confirming the transient nature of the steepening seen 
in Fig.[3]and Fig.|4] Ultimately, the density profiles will be approxi- 
mately unchanged from t heir initial state, in agree ment with a num- 
ber of other studies (e.g. lKazantzidis et al . 2004). In the righthand 
plot of Fig.[5] we display the evolution of potential energy with 
time for just the 5-20 kpc region. Notice how the potential well in- 
creases in depth after the interactions, confirming the steepening of 
the profile that is seen in Fig.[3] 

For both of the interaction scenarios displayed in Fig.[3]and 
Fig.|4] we also investigated the effect of varying the initial mass-to- 
light ratio of the satellite galaxy. In addition to the cases shown for 
an M/L = 10:1, we performed simulations using M/L = 20:1 and 
5: 1. In all cases there was very little variation in the results, and the 
overall effect of steepening followed by a return to isothermality 
was observed. We also examined higher mass interactions between 



Initial density profile 
Top solid line: Power-law fit t\ - 2.01 
Density profile after 3rd interaction 
Middle solid line: Power-law fit t\ = 2.70 
Density profile after 4th interaction 
Bottom solid line: Power-law fit t\ - 2.01 
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Figure 6. Density profile before and after a 3rd and 4th interaction between 
two 10 12 Mq objects. Here we have replaced the static potential (as in 
Fig.[3] and Fig. [4} with a particle based object. As in the static potential 
case, we continue to see the distinctive steepening followed by a return to 
an approximate isothermal state. 



10 12 and 10 13 Mq objects and again observed a similar pattern 
of steepening followed by a return to isothermality, although the 
larger of these two objects represents a mass normally only asso- 
ciated with central locations in clusters rather than in groups. To 
confirm that the exclusion of dynamical friction does not overly ef- 
fect the fundamental result, we replace the central static potential 
with a particle modelled object. We perform a simulation between 
two 10 12 Mq galaxies, i.e. the case in which dynamical friction 
should be at its greatest. As shown in Fig.[6]we continue to see the 
transient steepening in the density profile. 



4 DISCUSSION AND CONCLUSIONS 

We appear to require steeper than isothermal density profiles in 
order to reconcile certain lens galaxies with current estimates of 
the Hubble constant from HSTAVMAP. These galaxies include 
PG1 1 15+080, B 1600+434, HE2149-2745, and SBS 1520+530. One 
possible explanation is that certain early-type lens galaxies are off- 
set from the Fundamental Plane and as such early-type galaxies in 
general might display a heterogeneity in their structure and hence 
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mass density profiles. This theory proposes that the heterogene- 
ity results from group interactions between the lens galaxy and its 
fellow group members, resulting in central galaxies having isother- 
mal profiles and satellites having steeper than isothermal profiles. 
A number of studies seem to agree that the discrepancies in de- 
rived Ho values could be the result o f structural differences in 
lens galaxies as oppos ed to other factors Jfreu & Koopmansl20 02b ; 
iKochanek et~ai]|2006h . 

We have investigated group interactions between group mem- 
bers with masses ranging from 10 11 to 10 13 Mq. Our investiga- 
tions show a significant steepening of the density profile in the in- 
ner region, pertinent to strong lensing. This effect appears to be 
independent of the initial mass-to-light ratio, while the precise ini- 
tial orbital parameters made little difference to the fundamental re- 
sult provided consecutive interactions with the central galaxy took 
place. The magnitude of the steepening itself varied with the rel- 
ative masses of the interacting objects; for cases where the mass 
differed by an order of magnitude, as for the cases in interactions 
involving 10 11 vs. 10 12 Mq and 10 12 vs. 10 13 Mq, we observe 
originally isothermal profiles (rj = 2) steepen up to r\ = 3. For the 
cases where masses were of similar magnitude, i.e. an interaction 
between two 10 12 Mq group members, we observe a steepening 
of the profile slope up to ~ r\ = 2.7. In both cases, modelling such 
systems with isothermal models would result in a significant under- 
estimate of the Hubble constant. 

We stressed in §1 that most early-type galaxies appear to 
reside in groups jTonrv & Kochanekl ll999; Fassnach t et alj|2002l ; 
Willi ams et aT]|2006h . and as discussed bv lKochanek et all J2006t) . 
there is a play-off in lensing optical depth between satellite galaxies 
being more numerous than the most massive central group galax- 
ies, yet having a lower lensing cross-section (see also Oguri 2006). 
Even for satellite galaxies, the steepening in the inner region seems 
transient in nature, with consecutive interactions often returning the 
profile to an isothermal state within a timeframe of < 2 Gyr. This 
factor may help explain why lens galaxies that produce lower val- 
ues of Ho (i.e. those with possibly steeper profiles) are far fewer 
in number than those which agree with both the HST key project 
value for Ho and isothermality. In addition, it is important to note 
that this transient fluctuation in steepness is consistent with stud- 
ies showing that density profiles over longe r time scales can be ro- 
bust despite strong interactions or merg ers fkazantzidi s et al.l2004l ; 
iKazantzidis. Zentner & Kr avtsov 200g). 

We should highlight that other factors besides the environ- 
ment discussed in §1 are thought to contribute to discrepancies 
in derived Hubble cons tant values from certain systems. Recently, 
ISaha & Williams! feOOfih showed in an investigation with models 
of 35 galaxy lenses that shape-modelling degeneracies (e.g. caused 
by triaxiality) can also contribute to changes in the time delays and 
hence the derived Hubble constant value. This effect alone would 
need to be very significant to fully account for the systems with 
the lowest derived values of Ho- It seems probable that a combi- 
nation of effects, along with neglect of physical steepening due to 
stripping, can account for low Ho values. 

Looking towards the future, space and ground based obser- 
vatories such as GAIA and LSST will uncover roughly fifty times 
more multiply imaged quasars than we have identified today. With 
a greater sample of lens systems to study, and by combining lens 
models with stellar dynamical constraints, heterogeneity in early- 
type galaxies as a function of environment and redshift will be ex- 
plored in detail. This will also allow unprecedented constraints on 
the assembly and distribution of dark matter on galaxy and galaxy 
group scales. 
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